Maximal violation of the 13322 inequality using infinite dimensional quantum systems 
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The /3322 inequality is the simplest bipartite two-outcome Bell inequality beyond the Clauser- 
Horne-Shimony-Holt (CHSH) inequality, consisting of three two-outcome measurements per party. 
In ccise of the CHSH inequality the maximal quantum violation can already be attained with local 
two-dimensional quantum systems, however, there is no such evidence for the 13322 inequality. In 
this paper a family of measurement operators and states is given which enables us to attain the 
largest possible quantum value in an infinite dimensional Hilbert space. Further, it is conjectured 
that our construction is optimal in the sense that measuring finite dimensional quantum systems is 
not enough to achieve the true quantum maximum. We also describe an efficient iterative algorithm 
for computing quantum maximum of an arbitrary two-outcome Bell inequality in any given Hilbert 
space dimension. This algorithm played a key role to obtain our results for the 73322 inequality, and 
we also applied it to improve on our previous results concerning the maximum quantum violation 
of several bipartite two-outcome Bell inequalities with up to five settings per party. 

PACS numbers: 03.65.Ud, 03.67.-a 
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I. INTRODUCTION 



One of the most puzzling features of quantum theory is 
its nonlocal nature. Separated observers on a shared en- 
tangled state may carry out measurements on such a way 
that the correlations they generate are outside the set of 
common cause correlations In particular, such quan- 
tum correlations may find application in novel device- 
independent information tasks, which have no counter- 
parts in the classical world. They enable perfect security 
Q, randomness generation Q and state tomography Q 
without the need to trust the internal working of the de- 
vices. 

The concept of Bell inequalities is a particularly use- 
ful tool to detect nonlocal quantum correlations, as vio- 
lation of a single Bell inequality conclusively proves the 
nonlocal character of correlations. The standard scenario 
for a bipartite two-outcome Bell test is as follows. Two 
spacelike separated parties, Alice and Bob, both share 
copies of a quantum state jV') of a given dimension nxn. 
Alice (Bob) may choose between (ms) alternative 
measurements at random, where each measurement has 
two possible outcomes {0, 1}. In a single run of the ex- 
periment the correlations between the two {0, l}-valued 
observables Ai and Bj can be represented by the product 
AiBj. In order to obtain an accurate estimation of the 
correlations for each pair Alice and Bob repeat the 

experiment many times using a copy of the state lip) in 
each round. Averaging over many runs of the experiment 
yields the mean value (AiBj). 

The CHSH inequality Q is probably the most well- 
known and simplest example of a Bell inequality consist- 



ing of two measurement settings (m^ = = 2) both 
on Alice's and on Bob's part. Any correlations in the 
framework of local classical theories obey the following 
CHSH inequality [|, 

{AiBi) + {A1B2) + {A2B1) - {A2B2) - (Ai) - (Bi) < 0. 

(1) 

Curiously, quantum theory allows for a violation of the 
CHSH inequality, but the strength of nonlocal correla- 
tions is still limited, it obeys 0] 



(AiSi) + {A1B2) 
< I/V2- 1/2. 



{A2B1) - {A2B2) - {Ai) - (Bi) 



(2) 
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where now the expectation values can be expressed by 
{A,B,) =^{^\A, ® B,\^), (A,) = {i;\A, ® IbHj), (B,) = 
{iP\Ia ® i,j = 1,2 with {0, 1}- valued observables 

Ai, Bj. According to Tsirclson's theorem Q, the bound 
applies to any quantum correlations without making as- 
sumptions on the sort of measurements or the dimension- 
ality of the states involved. Though, the maximum value 
of I/V2- 1/2 can already be achieved with a maximally 
entangled two-qubit state. 

Inequalities, whose bound are saturated with quantum 
correlations without relying on dimensionality such as 
in Eq. ([2|) were coined as quantum Bell inequalities 
There exist various methods in the literature [8l-[l0|. to 
cite just a few, which enable one to derive quantum Bell 
inequalities. Also, several explicit constructions exist in 
the literature (e.g., [lTI - [l3 |). including so-called irrelevant 
ones [m (where the classical and quantum limits coin- 
cide). The method invented by Navascues, Pironio and 
Acm (NPA) [l3l is based on the solution of a hierarchy 
of scmidefinite programming (SDP) relaxations and is 
particularly useful since it gives better and better upper 
bounds on the maximum violation of an arbitrary Bell 
inequality by stepping to higher levels in the hierarchy. 
Moreover, the series of upper bounds in the hierarchy 
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keep to the exact quantum maximum in terms of com- 
muting measurements [iTI . [T6j . On the other hand, one 
can use heuristic algorithms to obtain nontrivial lower 
bounds in some finite dimensional Hilbcrt spaces on Bell 
inequalities, recovering the explicit form of the states and 
measurement operators as well. If the above computed 
upper and lower bounds coincide within numerical accu- 
racy for a given Bell inequality, then we may say that a 
quantum Bell inequality has been obtained, which delim- 
its the boundary of the quantum domain. 

In particular, in our previous papers [l3l [l7| we com- 
puted lower bounds on the maximum quantum viola- 
tions and determined the corresponding measurement 
operators and state vectors with two different methods, 
both involved some parametrization of the operators and 
applying a downhill simplex method to find the opti- 
mum parameters. The disadvantage of these methods 
is that the size of the Hilbert space was very limited. 
We could actually handle systems of maximum eight di- 
mensional real or six dimensional complex component 
spaces. Different sizes of component spaces require dif- 
ferent parametrizations for the operators, hence in order 
to extend these methods to a higher dimension would in- 
volve working out an appropriate parametrization. More- 
over, the choice of measurement operators were also lim- 
ited. For example, in the case of eight dimensional com- 
ponent spaces only projection operators projecting onto 
four dimensional subspaces were allowed. 

In spite of the limitations we could get the maximum 
quantum violation of the great majority of inequalities we 
considered (the list comprises 241 bipartite Bell inequal- 
ities with up to five settings per party collected from 
Refs. [13, [3 and detailed results concerning their 
optimum violations are presented in the web page [20|). 
However, there were still a few exceptions, where the up- 
per bound value resulting from the NPA method [l^l did 
not match the best lower bound result. The most in- 
teresting one was the case of 13322 ■ This is the smallest 
case we considered, and perhaps the simplest tight Bell 
inequality after the CHSH one, with only three measure- 
ment settings per party. It was introduced by Froissart 
[2l| back in 1981, and recently reinvented in Refs. [22l . [23j . 
It reads 

/3322 =-(^2) - (Bl) -2(^2) 

+ (AiBi) + {A1B2) + {A2B1) + {A2B2) 

- {A1B3) + (A2B3) - (^3^1) + (A3B2) < 0. (3) 

In a local classical model we have the maximum value of 
0, while the largest violation one could get with qubits 
was 0.25, which could already be achieved with a max- 
imally entangled pair of qubits (see e.g., [ill O, lisl - 
[isl . [23 )■ On the other hand, the best upper bounds 
are based on the NPA method [l^l and at level three it 
yields the significantly higher upper bound, 0.250 875 56 
We could even go above level three to an inter- 
mediate level in [isj , and presently we have got the upper 
bound 0.250 875 38 at level four. From the dependence 



of the bound on the level it was derived it seemed to be 
clear that it would not go much lower (note also that the 
computational complexity of the SDP problem increases 
dramatically with higher levels of relaxations). In fact, 
it is the above observation, which motivated us to search 
for a quantum violation beyond the two-qubit value of 
0.25. 

In the present paper we set out to resolve the puzzling 
problem concerning the maximum quantum violation of 
the /3322 inequality. To this end, we introduce in Sec. |TI] 
an efficient iterative algorithm which was applied to ex- 
plore the largest quantum violation of 13322 for states 
up to local dimension 20. Sec. IIIII contains the main re- 
sults, presenting the explicit construction of states and 
measurement operators, which give by means of an iter- 
ative method the conjectured maximal quantum value of 
-^3322 in function of the dimensionality. In particular, for 
dimensions very large we recover the upper bound com- 
puted with the NPA method at level four within high 
numerical accuracy, thereby establishing a quantum Bell 
inequality for 73322 ■ We note that the optimal quantum 
state in the infinite dimensional space is far from the 
maximally entangled state, thereby supporting the claim 
that entanglement and nonlocality are different resources 
[2^ . In Sec. IIVI wc investigate the remaining 19 Bell in- 
equalities from a set of 241 inequalities (plus two addi- 
tional symmetric 4-setting inequalities from Ref. [25|). 
where the tightness of the quantum bound could not be 
proven previously, and in some of the remaining cases 
now we manage to close the gap. Sec. fVl summarizes the 
results achieved. 



II. THE ITERATIVE ALGORITHM 

As mentioned above, the drawback of our previous 
heuristic methods in Refs. [l^ [l3l for computing lower 
bounds on the maximum quantum violation of Bell in- 
equalities is that the size of the local Hilbert space we 
could handle is limited up to dimension eight. With 
the iterative method we introduce here for two-outcome 
Bell inequalities wc can go to higher dimensions, which 
is only limited by the computational difficulties due to 
the increasing complexity of the problem with increas- 
ing Hilbert space dimensions. Similarly to the previous 
methods, this algorithm does not guarantee to converge 
to the global optimum solution. In case of larger spaces 
wc may miss it even after tens of thousands of restarts 
with different initial values. Nevertheless, it managed to 
find every optima we derived with our previous methods, 
and in almost all cases it has done it considerably faster. 

The problem to be solved consists of maximizing the 
quantum value of the Bell expression with coefficients 
M^i,, which can be written as: 

niA rriB 

Q = max^^AV(i^l4®B,|^), (4) 
^=0 1/=0 
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where there are niA and ms measurement settings for 
Alice and Bob, respectively, {1 < jj, < ttia) and B^, 
(l < ly < are the measurement operators of Alice 

and Bob, respectively, Aq = Ia and Bq = Is are unit 
operators in the component spaces, and \ip) is the state 
vector. The maximum can always be reached by a pure 
state. In case of two-outcome measurements it is enough 
to consider projection operators as the measurement op- 
erators [13] ■ When there are only two parties, a further 
simplification is that when we take a matrix representa- 
tion of the expression above, we may choose the bases 
in the component spaces such that they correspond to 
the Schmidt decomposition of the state vector j-f/;). If we 
confine ourselves to n-dimensional component spaces we 
get: 



Q„ = max ^ 5] ^ ^ M,^A>;^B-^^KX, 



(5) 



fj,=0 1^=0 i=l j=l 



where A^^ and B^^ arc components of matrices of A^ and 

Bi,, respectively, and A.^ are the Schmidt coefficients of 
the state vector. If n is smaller than the dimensional- 
ity of the component spaces required for the maximum 
quantum violation, then Q„ < Q. 

To solve the problem first we choose appropriate ran- 
dom matrices and numbers for Bob's measurement op- 
erators, and for Aj. To get the matrices, first we take a 
diagonal matrix with zero and one diagonal values chosen 
randomly, then we apply many random two dimensional 
unitary (or if we confine ourselves to real matrices, or- 
thogonal) transformations. For Ai we take positive num- 
bers between zero and one from a uniform distribution, 
then wc normalize them. In the iterative algorithm the 
first step is to calculate the optimal measurement opera- 
tors of Ahce, given Bob's operators and the state vector. 
This can be done directly. Equation ([5|) can be rewritten 
as: 



Q„ = max ^ ^ ^ A';^X^ = max ^ Tr{A,X,), 

fl—0 2—1 j—1 fl—0 



where 



(6) 
(7) 



i/=0 



We can get the matrix of the optimum (1 < M ^ itt-a) 
the following way. First we diagonalize the matrix of . 
Then we create a diagonal matrix. We choose its diag- 
onal matrix element one where the diagonalized matrix 
of Xf^ contains a positive number, and zero otherwise. 
Then we transform this matrix with the inverse of the 
transformation that diagonalized the matrix of X^. To 
show that this is the matrix of the optimal A^^ we note 
that trace is invariant to basis transformation, and it is 
easy to see that the optimum matrix of projector A^ in 
the basis diagonalizing is the one we have chosen. 



The second step is to derive Bob's optimal matrices 
while Alice's matrices and the state vector are fixed. This 
can be done equivalently to the first step. These two 
steps are the same as the steps of the see-saw algorithm 
by Werner and Wolf [2^, who used it to get maximum 
violation with fixed (not necessarily pure) states. A vari- 
ant of the see-saw algorithm for Bell inequalities with 
multiple outcomes was also devised [13, [l^ . 

In the third step the best state vector is calculated for 
the measurement operators. It is done the same way as 
in our previous method |13j . An algorithm applying this 
very third step has been also used by Liang et al. [1^ 
to compute quantum optima for multiple-outcome Bell 
expressions. From Eq. ([5]) it is clear that the optimal Ai 
corresponds to the eigenvector belonging to the largest 
eigenvalue of the matrix 



EE 



(8) 



and the eigenvalue is just the value of the Bell expres- 
sion. It is not sure that the components of the largest 
eigenvector will all be positive real numbers. If they are 
not, they can not be called Schmidt coefficients. How- 
ever, this does not affect anything, it is best just to leave 
them as they are. If we want to see the final result for 
the state vector represented by its Schmidt coefficients, 
it is enough to make the appropriate transformation at 
the end of the calculation. We simply have to multiply 
rows and columns of either Alice's or Bob's matrices with 
phase factors. 

The algorithm consists of repeating these three steps 
until the value of the Bell expression reach convergency. 
As this number may only increase in each step, conver- 
gence is ensured. What is not ensured is that we will 
converge to the global optimum of the problem. There- 
fore, we have to repeat the full procedure many times 
with different initial values. We may save computation 
time if we stop a run when it clearly goes towards an 
inferior solution. For example, if we are interested in a 
solution belonging to a larger Hilbert space, we may stop 
whenever too many components of the state vector falls 
below some threshold. 

It is not difficult to extend the method to more than 
two participants. In that case there is no Schmidt de- 
composition, which would ensure that one can choose 
the bases in the component spaces such that the n^- 
dimensional state vector has only n nonzero components. 
However, one can still use Eq. (jS]) to calculate the op- 
timum measurement operators for Alice, and analogous 
equations for the operators of the other participants. The 
only difference is that the formula for the matrix of Xf^ 
(and the corresponding matrices for the other partici- 
pants) will be somewhat more complicated than Eq. ([7|) , 
involving summations to more indices. The optimum 
state vector can also be calculated as the eigenvector be- 
longing to the largest eigenvalue of a matrix, but this 
will be a more general matrix in the full ??^-dimensional 
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Hilbcrt space. As the third step of the algorithm in this 
case involves the diagonalization of a much larger matrix 
than the first two steps, it is more efficient to do it less 
frequently, that is to repeat the first and the second steps 
several times before updating the state vector. 



III. FAMILY OF CONSTRUCTIONS 

From the optimality of formula ^ it follows that the 
matrices of the optimum Ai, A2, A^, Bi, B2 and are 
such that they satisfy the following Eq. dH), Eq. (fTO| . 
Eq. (HH), Eq. Eq. ^ and Eq. respectively: 



^Iji^lj + Bfj - Bf^)X^Xj 



^ A^^iBl^ + Bf^ + Bf^ - S,,)X,X, 



J2 4(4 + 4-4-<5.,)A,A, 



J2 s5(4+4+4-2<5,,)A.A, 



•t,i=l 



^4(4-4)A.A, 



max (9) 

max (10) 

max (11) 

max (12) 

max (13) 

max (14) 



Meanwhile, Ai are the components of the eigenvector be- 
longing to the maximum eigenvalue of the matrix 



M,;. 



(A 



2B 



(A 



Bl 



By) + Bf^{A^^~A^^). 



(15) 



The matrices of the measurement operators we got with 
the iterative method described in Section [IT] as the best 
solutions with all A^ differing from zero had very spe- 
cial forms. We found most of the matrix elements zero 
within expected numerical accuracy, and for odd number 
of dimensions it turned out that with some appropriate 
reordering of the bases we could also make all matrices 
blockdiagonal with one block having a single element, 
while the rest consisting of two by two blocks, with all 
matrix elements real. In the cases of Ai, A2 and -B3, the 
block of size one is the first one, while in the cases of 
Bl , B2 and A3 it is the last one. The actual form of the 



matrix of A2 may be written as 
/I 



A2 



1-C2 

2 

£2. 
2 



£2. 



2 

£i 
2 



£4 
2 

I+C4 
2 



V 



2 

Sn-l 

2 



1+cl 



(16) 

with Si = -^/l — c| arc all positive. The matrix of Ai 
tuns out to be the same, but the offdiagonal elements 
have a negative sign. The form of B2 is the following: 



Bo = 



l+ei 
2 

SI 

2 



SI 



I+C3 
2 

£3 
2 



£3 



Sn-2 
, 2 



v 



1 + C 

(17) 

The matrix element _B^j„ may take the value of 1 and 0, 
when c„ = and c„ = —1, respectively. For smaller di- 
mensions we got the former value for all best solutions, 
but it turned out that from dimensions larger than n = 79 
the other value gives the better solution. The notation 
of = 1 -I- Cn will become clear later. The matrix 
of Bl is the same as that of -B2, but = 0, and like 
in the case of Ai, the offdiagonal matrix elements are 
negative. Matrices A3 and B3 have even simpler forms: 

= Bf-^ = 1, while all elements of the two by two 
blocks are 1/2. It is easy to check that if all A; > 0, 
then the matrices of A3 and B3 do satisfy the optimality 
conditions Eq. ([TT|) and Eq. p4)). respectively. The ma- 
trix of A3 has the same block structure as B2 — Bi, as 
it should, and the optimality condition can be checked 
block by block. The value of the single element block 
in the matrix of B2 — Bi, that is — B^^^ is either 2 
or 0, depending on c„, in the first case the optimum is 

= 1, while in the other case the chosen value 1 is 
just as good as value would have been. The two by 
two blocks in the matrix of B^^ — have zero diago- 
nal elements, and positive nondiagonal elements. As all 
Ai have the same sign, the number they will be multi- 
plied with is also positive. Therefore, the optimum two 
by two block of the matrix of A3 is the symmetric two- 
dimensional matrix with eigenvalues and 1 having the 
largest possible nondiagonal elements. This matrix is the 
one whose all four elements are 1/2. 

From Eq. pU)) we can get equations determining the 
parameters of the matrix of A2 that is Ci with i even (see 
Eq. (dni) in terms of the parameters of the matrices of 
Bob's measurement operators, that is c; with i odd. and 



+ Cn J 



the Schmidt coefficients A,. As the matrix oi Bi + 132 — I 
is diagonal, the block structure of Ai is the same as that 
of B3. The matrix oi Bi + B2 + B3 - I is; 

Si + B2 + S3 - / = 



h- ci 
^ 1 
2 



C3 



C3 



1 

1 2 
5 + 



1 — ^ 1 

2 '-"-2 2 



The value of is determined by the sign of ci, the 
A\i = 1 value is correct if ci > 0, which turns out to 
be true. Each Ci, where i is even, occurs only in one 
two by two block of the matrix of A2 (see Eq. 
Its optimum value depends on the corresponding block 
in Eq. (fTSi) . containing d-i and Ci+i. We denoted i?,^j„ 
as 1 + c„ such that the last block has the same form 
as the others. Using Eq. ([TU|) we arrive at the equation 
determining the optimum value of c^: 



Ct-i A: 



(1 



1 



+ Ci+i Xi+i 



(19) 



Multiplying the equation by four, dropping terms con- 
stant in Ci and substituting Si = y^l — cf we get: 



c,[(l + 2c,+i)A2^i - (1 - 2c,_i)A2] + 2^l-c^\\+i 
^ max, (20) 

from which it follows that 

(21) 



where 



^+4A^A^i 



r, s (1 + 2q+i)A2^i - (1 - 2q_i)A2. (22) 



With Ci chosen this way Eq. the optimality condition 
for Ai is also satisfied. If we started from this condition 
instead of the one for A2 , we would have arrived at exacly 
the same results. 



In the same way as above, from Eq. and Eq. 
we can derive formulae for the optimum values of the 
parameters of the matrices of Bi and B2, that is for Ci, 
where i is odd. These formulae turn out to have exactly 
the same form as the ones derived for even i values, that 
is Eqs. (|19ll22p are valid for odd i. To get ci, we have 



(18) 
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FIG. 1: Distance of Q„ from tlie upper bound of Qmax ~ 
0.250 875 38 is shown as function of the dimensionality for 
cases of c„ = and c„ = —1. 



to introduce the notation cq = 1. The value B^^ = is 
the optimum value. Either — (that is c„ = — 1) or 
^nn = 1 (that is c„ = 0) is a possible choice. Consistency 
with the optimality condition is fulfilled if c„_i < in 
the first case, and c„_i > in the second one. 

With the specific forms of the measurement operators 
above, the matrix M of Eq. whose eigenvector cor- 
responding to its largest eigenvalue gives the Schmidt 
coefficients Ai, is tridiagonal, and its nonzero matrix ele- 
ments are: 



Mi, 



Ci-lCi + 



Cl-l 



M, 



(«+i) 



2 



(23) 



We may apply the iterative method with confining our- 
selves to the present specific forms of the operators, and 
get solutions for very large dimensions. This way we may 
increase the dimensionality of the component Hilbert 
spaces even to a thousand, or more. The distance of the 
maximum value of the Bell expression from the upper 
bound 0.250 875 38 calculated at level four is shown in 
Fig. [T] as a function of the dimensionality, both for c„ = 
and c„ = —1. For lower dimensions the former family of 
solutions is much better. Still, even this gives a value 
larger than 0.25, which one can obtain with a pair of real 
qubits only if n > 12. Unfortunately, this family does 
not converge to the upper bound. The other family with 
c„ = — 1 is initially inferior, but above n = 79 it over- 
takes the other one, and converges to 0.250 875 384 514, a 
value well consistent with the upper bound: it is slightly 
larger, but within the numerical accuracy of the bound. 
In Fig. [5] and Fig. [3] we show Ai and Ci for n = 99 with 
c„ = and c„ = —1, respectively. For small i the Ci 
is close to co=0, and tends to a constant value fast. In 
the case of c„ = it stays basically constant up to the 
last few point, when it drops towards the zero value of c„ 
fast. When c„ = — 1 and n > 19 there are more than one 
solutions. In all solutions c; > for small i and Ci < for 
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FIG. 2: The values of A; and c; are shown for 71 = 99 with 
c„ = 0. 



FIG. 3: The values of \i and Ci are shown for n = 99 with 
Cn = —1. 



large i. In the different solutions the change of sign hap- 
pens between two neighbouring integers. For very large 
n it always happens above 16 and below n — 16. Among 
all these solutions the best is the one in which Cj changes 
sign exactly in the middle, that is between (n— l)/2 and 
(n + l)/2. In Fig. Owe show the result for this solution. 
In particular, the upper curve in Fig. [3] shows the be- 
haviour of Ci. The curve starts like in the case of c„ = 0, 
that is it drops from cq = 1 to a constant value fast, that 
stays there until it goes to minus one times the constant 
value in a relatively short interval, and then it stays al- 
most constant again until it turns towards c„ ~ —1 at 
the end. Where is constant, behaves exponentially. 



make the change of size to happen at the middle, that is 
we enforce c„/2 = Oi for large n the result will converge 
to the same suboptimal value as in the case of the c„ = 
family. Probably the reason this family fails to converge 
to the true optimum is that c, must approach zero at an 
integer value, instead of midway between two integers. 

In the matrices of the measurement operators we made 
very specific choices for the values of the single element 
blocks, leaving only c„ as a single free parameter. We 
note that there are other choices consistent with the op- 
timality conditions, but they lead to solutions equivalent 
with the ones discussed. 



Knowing the behaviour of c,;, it is easy to choose initial 
values such that the iterative procedure will converge to 
the desired solution for the first time, there is no need 
for repeated runs. For example, we may initialize Ci as 
+C for i < n/2 and -C for i > n/2, with C = 0.9. Then 
we can start with the third step instead of the first one, 
that is with the determination of A^ . 



The families of solutions may be extended to even val- 
ues of n. The matrices of the measurement operators are 
blockdiagonal, like for odd n. The matrices of Bi, B2 and 
A3 consist of n/2 two by two blocks. The blocks of Ai, 
A2 and are also two-dimensional except for the first 
and the last ones, which have single elements. The actual 
form of the matrices can be parametrized with Ci anal- 
ogously to the odd n case, which will satisfy the same 
equations as optimality conditions, and the eigenvalue 
equation for will also have the same form. Figure [T] 
does contain the results for n even. The curve for Ci in 
the case of the best solution with c„ = — 1 is not fully 
symmetric. It changes cither very nearly half way be- 
tween n/2 — 1 and n/2 or between n/2 and n/2 -I- 1 . If we 



IV. OTHER EXAMPLES 

In Ref. [l3| we calculated both lower and upper bounds 
for all known tight binary bipartite Bell inequalities with 
up to five measurement settings per party. For most cases 
we managed to find explicit solutions saturating the up- 
per bound. There were 20 exceptions. The case of 13322 
was discussed in the previous chapter. Now we have ap- 
plied the iterative algorithm with up to 25 dimensional 
component spaces for the remaining 19 problems. In 11 
cases we have found better solution than the ones we re- 
ported in Ref. [l^ . The solution for Ags could have been 
found by our previous methods, but we missed it. All the 
other ones involve larger Hilbert spaces than those meth- 
ods could handle. We have also made the upper bound 
tighter by going to higher levels by using a computer 
with larger memory. As previously, the calculations have 
been done with Borcher's code CSDP for semidcfinite 
programming (36| . For cases with five measurement set- 
tings per party we could do the calculation on level two 
plus aa'b plus abb'. For the smaller cases we could afford 
to go up to level three, except for 13322, where we did level 
four. The notion of levels and partial levels, and their no- 
tation is explained in Ref. [ll| , and also in Ref [l^ . This 
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TABLE I: Results for the 19 Bell inequalities for which the tightness of the quantum bound has not been proven. The 
local Hilbert space necessary to achieve the best lower bound is denoted by H, column 'New' tells if this lower bound is an 
improvement over our previous result, and column 'Level' is the level to derive the upper bound. 



Case Lower bound H New Upper bound- Level 

lower bound 



7-20 
-'4422 


0, 


.4676794 




no 


0, 


.0000000 


L3 


+ 


aa'bb' 




t30 
-'4422 


0, 


.4363842 


R" 


yes 


0, 


.0008787 


L3 








i34 
-'4422 


0, 


.4492657 


R* 


no 


0, 


.0000000 


L3 


+ 


aa'bb' 




i48 
-'4422 


0, 


.7516220 


R^2 


yes 


0, 


.0000120 


i3 


+ 


aa'bb' 




i62 
-'4422 


0, 


.7500755 


R« 


no 


0, 


.0000166 


L3 








i73 
-'4422 


1, 


.1223170 


r20 


yes 


0, 


.0002409 


L3 








7-97 

-'4422 


1, 


.1583626 


R2 


no 


0, 


.0001425 


L3 








Ai4 


0, 


.4759513 


r20 


yes 


0, 


.0018953 


L3 








All 


0, 


.3260601 




yes 


0, 


.0001176 


L3 








A47 


0, 


.4608544 




no 


0, 


.0020847 


L2 


+ 


aa'b + 


abb' 


Ae2 


0, 


.4065268 


R^^ 


yes 


0, 


.0003457 


L2 


+ 


aa'b + 


abb' 


Ae4 


0, 


.3900890 


R^ 


no 


0, 


.0013878 


L2 


+ 


aa b + 


abb' 


Aes 


0, 


.3688996 


R^ 


yes 


0, 


.0000000 


L2 


+ 


aa'b + 


abb 


Ae,7 


0, 


.3990671 


R^ 


no 


0, 


.0000000 


L3 








Aes 


0, 


.4025522 


Ri« 


yes 


0, 


.0024807 


L2 


+ 


aa b + 


abb' 


Aso 


0, 


.3769863 


R* 


no 


0, 


.0002557 


L2 


+ 


aa'b + 


abb' 


As2 


0, 


.4708838 


r20 


yes 


0, 


.0000341 


L2 


+ 


aa'b + 


abb' 


As4 


0, 


.6352087 


Ri« 


yes 


0, 


.0000747 


L2 


+ 


aa b + 


abb' 


As9 


0, 


.3035637 


R" 


yes 


0, 


.0022575 


L2 


+ 


aa b + 


abb' 



way the two bounds met for Ae5. We sent Aer at level 
three, and J4I22 and J4422 at level three plus aa'bb' to 
Brian Borchers, as test cases for his new code which re- 
quires much less memory than CSDP for very large cases. 
For Aqj this new upper bound does agree with the lower 
one. Unfortunately, for J4422 there remained a difference 
of 0.000012, significantly less than the level three value 
of 0.0000657, but more than numerical uncertainty. At 
an even higher level it would probably disappear, but it 
is not sure. Therefore, we can not state that our solu- 
tion with 12-dimensional component spaces is already the 
optimum one. The case J4J22 at level three plus aa'bb' 
seems to be a difficult case for the new code, it could 
solve it only at a reduced accuracy, and the value it has 
given for the upper bound turned out to be lower than the 
lower bound. Therefore, very probably the solution with 
a pair of real qubits we found the best so far is actually 
the optimum. If we added to level three just one quarter 
of the aa'bb' type terms, we could use CSDP. In the cases 
of /I422 and J4I22 this was already enough to reach the 
lower bound. By looking at the dependence of the upper 
bound on the level it has been calculated, we are almost 
sure that for JI422, -^^4422! ^47, ^62, ^64, ^82 and Agg the 
present lower bound will not be reached, therefore we do 
not yet have the optimum solution. Probably 73322 is not 
the only inequality requiring infinite dimensional Hilbert 
spaces for maximum violation. Actually, the difference 
in the cases of J4I22 and Ag,2 is very small, but it docs 



not change much with increasing level. For Ago and As4 
the upper bound would more probably converge to our 
solution, while for the rest of the cases we can not tell 
from the present results. 

Bancal et al. j25j found recently two symmetric in- 
equalities with four settings per party, which they called 
S'51 and 552. They were also included in our calculations. 
In fact, both of them could be saturated with measure- 
ments acting on qubits, we found the maximum violation 
to be 1.0135274 and 0.87038004 for S'51 and S52, respec- 
tively. 

The details of our best solutions, including Schmidt co- 
efficients and the matrices of the measurement operators 
are on our web site [20ll. 



V. CONCLUSION 

To summarize, we investigated the 13322 Bell inequal- 
ity and its maximal quantum violation. Although 13322 
is one of the simplest Bell inequalities with just three 
measurements per party, surprisingly we found that the 
maximum value of 0.25 achievable with a pair of qubits 
could be only be overcome by using states of dimension 
at least 12 x 12. Moreover, we found a family of measure- 
ment operators and states, which attains numerically in 
very large dimensional Hilbert spaces (and hence in the 
limit of infinite dimensions as well) the largest possible 
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value of the Bell expression 13322 allowed by quantum 
theory. Taking as a conjecture that our construction is 
optimal (in the sense that no finite dimension suffices to 
achieve the true quantum maximum for 13322), to our 
knowledge this constitutes the first example of a Bell 
scenario concerning finite number of measurements and 
outcomes, where measuring finite dimensional quantum 
systems is not enough to obtain exactly all quantum cor- 
relations. In this respect we wish to invoke the concept 
of dimension witnesses introduced recently in Ref. [30|. 
In a broader sense, the aim of this concept is to put a 
lower bound on the Hilbert space dimension needed to 
reproduce the statistics arising from a Bell experiment. 
Dimension witnesses have been found or conjectured in 
many different scenarios, such as for the bipartite setting 
[il Q [H miiSai] for the multipartite setting [H and 
even for the case of a single system [35| . 

In particular, for measurements with binary outcomes 
it was shown how to get dimension witnesses analyti- 
cally for any dimensions [3l| . Actually, these results en- 
tail that no finite dimension is sufficient to generate the 
whole set of quantum correlations provided an arbitrary 
number of two-outcome measurements is involved. Con- 
sider now a scenario where the number of measurements 
(each having a finite number of outcomes) is finite. Does 
there exist a Bell inequality requiring infinitely large en- 
tangled states for obtaining the maximum possible quan- 
tum violation? That is, are there quantum correlations 
corresponding to this scheme attainable exactly by mea- 



suring infinite dimensional quantum systems? This is a 
question posed recently by Navascues et al. [ll| . 

Indeed, in the present paper we suggested, supported 
by thorough numerical computations, that no finite di- 
mension suffices to maximally violate the 13322 inequality. 
While we cannot guarantee that for any finite dimension 
there exist no other construction beyond ours giving a 
larger quantum value, nevertheless we showed that our 
best conjectured values tend to the upp er limit computed 
with the aid of the NPA method We pose it as a 

challenge to provide an analytical proof of our conjec- 
ture. Also, in our view it would be interesting to find 
other, bipartite or multipartite. Bell inequalities beyond 
^3322 (but still involving finite number of settings and 
outcomes) requiring infinitely large entangled states for 
their maximal quantum violation. 
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